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Abstract 



We present the Hit Backspace Once More (HBOM) method for 
correcting a measurement for the effect of an imperfect detector. The 
HBOM method is a model-independent and potentially data-driven 
technique that repeatedly applies a parameterisation of the detector 
effect to observed data. The correction is determined by extrapolating 
the data so-obtained to a detector effect of zero. We demonstrate 
this technique using the two particle correlation function, which is 
an observable that can otherwise be difficult to correct for systematic 
shifts introduced by the detector. 

1 Introduction 

A measurement made on data from a collider experiment will normally in- 
clude artefacts arising from the interaction of the final state particles with 
the detector. Such "detector effects" include energy and angular smearing, 
production of additional particles through secondary scatterings, and ineffi- 
ciencies in the detection of particles. A measurement that contains detector 
artefacts is of less use than one that does not; in order to correctly com- 
pare such a result to either theoretical predictions or other measurements of 
the same observable, one needs access to a full description of the detector 
effects. Detailed knowledge about the detector, typically in the form of a 
sophisticated simulation package, usually only exists within the experimen- 
tal collaboration responsible. There is no guarantee, or indeed expectation, 
that the detailed detector simulation will be available many years from now. 

A good measurement will, therefore, de-convolute the detector effects 
from the fundamental physics processes and publish results that have been 
unfolded to match what one would observe with a "perfect" detector. Such 



results are directly comparable to a hadron-level Monte Carlo (MC) simula- 
tion with no need for a detector simulation. 

Several methods for unfolding the detector effects exist. One of the sim- 
plest is bin-by-bin unfolding, for which the measurement is performed on 
an ensemble of Monte Carlo simulations both with and without a detector 
simulation. Correction factors are determined for each bin of the result from 
the ratio of the value with and without the detector simulation. Bin-by-bin 
unfolding has the advantages that it is easy to understand and can work 
quite effectively on simple observables. However, bin-by-bin unfolding does 
introduce a dependence on the Monte Carlo model (s) used to derive the 
bin-to-bin migration correction factors and, especially for more complicated 
observables, it tends to produce a bias in the shape of distributions. 

Bayesian unfolding [1] expresses the probability that, given an observed 
value Xf s , the true value is Xf ue . This expression requires knowledge of 
the set of probabilities that a true value _Xj™ e will be observed as Xf bs , 
together with a prior model for the probability distribution of the true val- 
ues. An iterative procedure is applied within Bayesian unfolding such that 
the corrected data from the first iteration is used as the prior for a second 
iteration. 

Iterative Bayesian unfolding has the advantages that it works even for 
complex observables and in principle, for a sufficiently large statistical sam- 
ple, it contains negligible model dependence. However, the iterative sequence 
may take a long time to converge, which can introduce large fluctuations into 
the result. The choice of prior can also have an effect on the unfolded distri- 
butions if there are an insufficient number of iterations. Bayesian unfolding 
is also quite an involved procedure and can be conceptually hard. This latter 
point is not insignificant because in order to properly assign to a measure- 
ment any systematic uncertainties arising from the unfolding procedure, one 
needs to fully understand the process. 

Any unfolding procedure is simply trying to quantify as well as possible 
the effect that the detector has on an observable. Conceptually, one may 
consider a detector to be a (very complicated) mathematical operator, A, 
that has been applied to the observable, Xf X° bs = X { (0) = AX\ rue = 
AXi ( — 1). If A, or a reasonable approximation to it, is known then it may 
be applied again to X { (0) to give X; (1) = AX { (0) = AAX\ rue . In this way 
(and without knowing X{ (—1)) a sequence of identical observables with an 
increasingly large detector effect can be constructed: Xi (0) , Xi (1) , Xi (2) ... 
X{ (N). Under the assumption that Xi (iV) is a smooth function of N, a 
numerical fit can be made to the sequence of Xi. Such a fit can be evaluated 
at N = — 1 to provide an estimate of the true value of the observable in the 
absence of any detector effects. 

In this note we show an example of this method applied to a two particle 
correlations observable measured in an ensemble of Monte Carlo generated 
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events. Detector effects can often express themselves in a complicated way for 
correlation observables because the loss of a single particle from the detector 
makes a contribution to the observable that depends on all of the remaining 
particles in the same collision event. As such, bin-by-bin unfolding does not 
work to correct the two particle correlation function because it introduces a 
significant bias to the shape of the function. 

2 Observable Definition 

The two particle correlation function has been used by several experiments 
in order to investigate soft non-perturbative effects and particle production 
mechanisms [21 El SI IS]- The definition used here was fully inclusive for 
all events containing more than two charged particles. All charged parti- 
cles whose transverse momentum and pseudo-rapidity {jpr and rj, respec- 
tively) satisfied pt > 100 MeV and \rj\ < 2.5 were accepted into the analysis. 
The correlation function was defined in terms of a foreground distribution, 
F (At), and a background distribution, B (At), where At was the separation 
between a pair of particles in either r\ (Ar\) or azimuth (A(f>). As such, the 
two particle correlation function probed the correlations between pairs of 
particles as a function of their At] or A0 separation. 

The foreground distribution was obtained by determining the At for each 
pair of charged particles within the same event and filling a histogram with 
those At values weighted by 2/N c h, where N c h is the multiplicity of charged 
particles in each event. The multiplicity-dependent weighting factor is a 
standard feature of two particle correlations that gives each event a weight 
according to (N c h — 1) /N p , where N p is the number of charged particle pairs 
per event. The factor of (N c h — 1) is present because the strength of corre- 
lations between particles was historically found to be approximately propor- 
tional to (N ch - 1) 0. 

The background distribution used was similar to the foreground distri- 
bution, but pairs of particles were taken from independent events before 
determining the At distribution. Unlike the foreground distribution, no 
multiplicity dependent weighting factor was used for the background. The 
background distribution was instead normalised to unit integral. By obtain- 
ing the correlations between particles in independent events, the background 
distribution reveals the contribution to correlations that arises solely from 
acceptance cuts (or, in a real experiment, effects caused by the detector). 

Having obtained the foreground and background distributions, the two 
particle correlation function, R(At), was produced as given by equation [T] 



where {N c h — l) evts is the value of (N c h — 1) averaged over the ensemble of 
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all events. The general shape of R(Arj) and R(A<p) is shown in figure [2] of 
section [5] 

3 Monte Carlo Sample 

A sample of 50 million minimum bias Monte Carlo events was generated 
with tune "4C" of the Pythia 8 Monte Carlo generator [6] . Tune 4C was pro- 
duced after a comparison of Pythia 8 to early results from the Large Hadron 
Collider, which motivated a somewhat reduced cross section for single and 
double diffraction as well as an improved description of the average trans- 
verse momentum per charged particle [7]. Mixed samples were produced 
containing non-diffractive, single and double-diffractive proton-proton colli- 
sion events at a centre of mass energy of 7 TeV. Events were analysed using 
the Rivet [8] Monte Carlo analysis framework. 

4 Pseudo-detector Simulation 

The effect of a detector on the observed charged particles was parameterised 
as a function of a particle's pr and 77. A detector will typically have an effi- 
ciency for detecting charged particles that rises rapidly from around 50 MeV 
to nearly unit efficiency at around 1 GeV or more. The efficiency for de- 
tecting a charged particle will typically fall gently as r\ departs from zero. 
The detector efficiency employed for this study was therefore in the form of 
an inverse tangent function of and a broad Gaussian in 77, as given in 
equation [2] 

*<™>-§~- , (™H-te)') < 2 > 

The parameters a p , pto anci ®r\ determine how quickly the efficiency rises and 
falls with pj- and r\ respectively. Realistic nominal values of a p = 0.15 GeV, 
Pto = 80 MeV and a v = 6 were used, which resulted in the efficiency shown in 
figure[T] In order to investigate the effect of using a detector parameterisation 
that does not perfectly match the true detector effect (if, for example, the 
particle detection efficiency were not perfectly known or were mis-modelled) 
a second (henceforth the alternative) detector parameterisation was defined, 
for which values of a p and were taken as 0.18 GeV and 7, respectively. 

The particle detection efficiency of equation [2] was applied to the Monte 
Carlo sample by generating a unique random number in the range {0, 1} for 
each particle in each event. If the random number for a given particle was 
greater than its detection efficiency then that particle was removed from the 
sample. This procedure was used for the initial generation of X, (0) from 
j£true an j £ or j. kg subsequent generation of the sequence X{ (1) ...Xi (N). 



4 




Nominal Efficiency 

Alternative Efficiency 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

1 



h- 

Q_ 

LU 1 




0.8 




0.6 


/.'" Nominal Efficiency 


0.4 


/.■' Alternative Efficiency 


0.2 





Figure 1: The particle detection 
efficiency as a function of pt and 
rj. Top left shows the 2D pt and 
rj dependence, while the top right 
and bottom left show the sepa- 
rate pt and r] dependence, respec- 
tively. The nominal and alter- 
native detector parameterisations 
are shown for the latter two plots. 
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In addition to employing an incorrect particle detection efficiency, the de- 
tector effect, A, that is repeatedly applied to Xj (0) to obtain Xj (1) ..Xj (N) 
may lack certain aspects of the detector's behaviour. In order to investigate 
this possibility, we also introduced additional particles into a sample of events 
to approximate the effect of secondary scatterings within the detector, but 
without including that effect in A. The probability that a particle splits into 
a pair of particles was parameterised as a function of rj as in equation [3] 

P(r ] )=e W ^y^j -0.97 (3) 

The resulting pair of particles was given an opening angle 9 S with a proba- 
bility P (Q s ) given in equation [4] 

P{0.) = yfl (l + tan 2 (|)) exp (-Stan 2 (|)) (4) 

In general, a detector simulation should also include smearing of the parti- 
cle energies and angles. However, in this case, the observable did not have 
a strong dependence on particle energies because only a charged particle's 
rj and 4> coordinate enter into the two-particle correlation function. Thus 
energy smearing would only have a small effect due to the px > 100 MeV 
charged particle cut, and as such was not considered here. In any real ex- 
periment the observable measured here would be obtained from tracking de- 
tectors, which typically provide excellent angular resolutions of better than 
0.1 milliradians. Angular smearing was therefore also not considered, since 
any effect would be small in a realistic scenario. 

5 Application of the HBOM Method. 

The nominal detector parameterisation was applied to a sample of hadron- 
level truth events, from which the two particle correlation observable was 
then calculated. The result is an approximation to what would be observed 
if the true events were detected in a real detector. The same detector pa- 
rameterisation was then applied a further five times to the same ensemble of 
events to give Xj (1) ...Xj (5). 

The set of observables Xj (0) ...Xj (5) (each expressing an increasingly 
large detector effect) is shown together with the true value of the observable 
in figure [2] Note how the structure in the observable (the peak at Arj = and 
the gentle rise at high Ar/) is flattened out with each additional application 
of the detector parameterisation. Flattening occurred because, in removing 
particles from events, the detector effect weakened the correlations observed. 
The amount of flattening decreases with each subsequent application of the 
detector. In order to arrive at an estimate of the true observable, such 
flattening had to be reversed. Conversely, the Acp dependence of the two 
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particle correlation function develops a somewhat sharper peak at Acf) = 
(and correspondingly broader peak at Acf) = if) as the detector effect is 
increased. The different behaviour shown by the Arj and Acf> dependence 
of the two particle correlation function illustrates why a simple re-scaling of 
the results tends not to work for complex observables. 
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Figure 2: The A?] (a) and Acf) (b) dependence of the two particle correlation 
observable. The true value of the observable is shown as solid black circles, 
while the value observed by the parameterised nominal detector is shown as 
purple open triangles. Further applications of the nominal detector to the 
observable show an increasing flattening of the correlation function in Ar] 
and a slight sharpening of the correlation in Acf) near zero. 



The value of the observable in each of its bins was plotted as a function 
of the number of times the detector parameterisation had been applied. A 
polynomial of degree four was fitted to the resulting "iterations curve" for 
each bin of the observable. Figure [3] shows the fits for three different points 
in the two particle correlation function at Arj = 0, Ar] =1.6 and Ar] = 4.5. 
The values of Ar] shown here were chosen in order to be representative of the 
whole range of correlation function values. Figure [3] also shows additional 
fits that are discussed in section The true hadron-level value of the corre- 
lation function is shown in figure [3] (but that point was not used in the fit), 
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together with the estimate of that true value obtained from the fit and its 
corresponding fit uncertainty. The agreement between the point estimated 
by the fit at TV = —1 and the true value of the observable was found to be 
generally good. 
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Figure 3: Polynomial fits to the value 
of the two particle correlation function 
Vs. the number of applications of the 
detector, N. The value measured by the 
detector is at N = 0, the true value is 
shown by a star at N = — 1. The ob- 
served N = value is estimated from 
both a cubic and quartic fit to the N = 1 
and higher points, while the true N = 
— 1 value is estimated from a quartic fit 
to the N = and higher points. Fits are 
shown for three A77 values of (a), 1.6 
(b) and 4.5 (c). 



6 Results of the Correction Procedure 

The true value of the two particle correlation function is shown in figure [4] 
together with the value observed by the nominal detector and the corrected 
nominal detector. The correction has removed the vast majority of the nom- 
inal detector effect, although a small residual detector effect remains. Note 
that while the correction is quite large for the Ar) dependence, the same 
procedure applied to the A<p dependence provides a much smaller correc- 
tion that acts in the opposite direction; enhancing the peak at Ar/ = 
while slightly flattening the peak at Acj) = 0. The bottom panels of figure 
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[4] show the difference between the true and corrected values of the two par- 
ticle correlation function. The correction to the nominal detector was also 
determined by using the alternative detector parameterisation to construct 
Xi (1) ...Xi (5) (whilst using the nominal parameterisation to construct the 
detector-level observable Xi (0)). The corrected value thus obtained by us- 
ing a mis-modelled detector to apply the correction is also shown in figure 
[4] In this example, when using a polynomial fit of degree four to the detec- 
tor iterations curve, the mis-modelled alternative detector in-fact shows a 
slightly better agreement between the corrected and true observables. This 
better agreement occurs because the alternative parameterisation represents 
a stronger detector effect; a particle has an increased probability of being 
removed from the sample when using the alternative as opposed to nomi- 
nal detector. The alternative correction therefore over-corrects relative to 
the nominal correction, however this effect cancels with the slight under- 
correction provided by the nominal correction when using this particular 
polynomial fit to the iterations curve. 

The additional effects of particle splitting due to possible interactions 
with the detector material (equations [3] and Q are shown in figure [5] Unlike 
the inefficiency for particle detection exhibited by the nominal detector, the 
addition of such secondary particles had a much stronger effect on the A<p 
dependence of the correlation function than on the Ar/ dependence. In both 
cases the addition of particle splitting enhances the peak of the correlation 
function. This effect is easily understood because the particles so-added are 
necessarily correlated. 

Applying the the splitting effect to the N = detector level observable 
while using only the nominal detector to derive the correction led to a clear 
bias, as is shown in figure [5] The main effect is that, after correction, the 
central peak at Arj = and A(fi = is noticeably too high when using the 
nominal correction with the sample containing split particles. This happens 
because the effect of splitting acts in the opposite direction to the effect of 
the loss of particles, therefore correcting only for the loss of particles (as 
is the case with the nominal correction) leads to a correlation that is too 
strong. 

The splitting was corrected for by including the splitting parameterisa- 
tion of equations [3] and [4] in the detector function used to derive the ./V = 1...5 
iterations. Fitting a polynomial to the resulting iterations curve and extrap- 
olating to N = —1 corrected for both the tracks missing from the sample 
and the tracks introduced by the splitting. The corrected two particle corre- 
lation function so obtained shows an agreement with the true value that is 
as good as the corrected nominal sample, thus demonstrating that the effect 
of splitting can be entirely eliminated with this method. 
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Figure 4: Corrected and uncorrected Arj (a) and Acft (b) dependence of 
the two particle correlation function. Corrections were applied using linear 
polynomial fits of degree four to the iterations curves (see figure [3| . The 
nominal detector-level observable is shown as filled circles (colour online). 
The nominal detector after corrections is shown as un-filled circles, while 
a correction that uses the miss-modelled alternative detector is shown as 
un-filled triangles. The true value is shown as un-filled squares and the 
differences between the true and corrected values are shown in the bottom 
panel. 



7 Variation of the Fits to the Detector Iterations 

Fitting a polynomial of degree four produced a satisfactory correction that 
quite closely matches the generated hadron-level distributions. In order to 
evaluate the success of the fit, other possible fitting functions were also ex- 
plored. 

As a test of the validity of the fitted function, polynomial fits of degree 
three and four were made to the points for which N > 1 and extrapolated 
to the N = observed point. Such a test would be possible during a real 
experiment, in which the true value of the observable would be unknown. 
The resulting fits are shown in figure [3] while the estimate of the N = 
observed uncorrected value of the correlation function is shown in figure [6] 
Figure [6] shows that, in this case, the fit using the polynomial of degree four 
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generally exhibited less bias than the polynomial of degree three. Note that 
while a polynomial of degree four was guaranteed to fit five data points more 
successfully than a polynomial of degree three, there was no guarantee that 
the extrapolation to N = — 1 would be closer to the observed value. 

The iterations curve was very often (but not always) described reasonably 
well by an exponential decay because the number of particles removed at each 
detector application was proportional to the number of particles remaining 
in the sample. Polynomial curves of degree three and four were therefore also 
fitted to the logarithm of the absolute value of each point in the iterations 
curve. Note that, for a small number of points where the correlation function 
was near zero, logarithmic fits were not possible because after applying the 
detector effect several times, the correlation function can change sign. When 
plotted in this way on a log scale where possible, the value of the observable 
was often close to linear in the number of detector applications. The test of 
estimating the observed N = point from the N > 1 points was performed 
and is shown in figure [7] for the same Ar] values as figure [3j Fitting to the 
logarithm of the data points showed a generally good agreement between the 
estimated and true values, which was quite often better than the linear fit. 
However, fitting to the logarithm of the data could, in some cases, lead to 
an unstable result due to large contributions from the higher order terms of 
the polynomial. For such points it is better to use a fit without taking the 
logarithm of the data. 

The corrected value of the correlation function using a logarithmic fit 
is shown in figure |8j The logarithmic fit resulted in corrected curves that 
exhibit somewhat less bias than the linear fits of figure |4j although the loga- 
rithmic fits tend to show slightly larger uncertainties on the fit parameters. 

8 Discussion and Summary 

The HBOM method for correcting measurements for detector effects has the 
following advantages compared to other methods: 

• It is data-driven and model independent. A reasonable approximation 
to the detector behaviour is sufficient. 

• It can be generalised to a large number of observables, including multi- 
dimensional observables. 

• Complicated observables, including those affected by correlations, may 
be corrected. 

• Histogram bin widths do not need to be optimised to avoid migrations. 

• No prior physics model is required. 

• The method is conceptually simple. 
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The method is ultimately limited by the approximation to and parame- 
terisation of the detector response that is applied to the data. In principle, 
a model of the detector may be constructed without reference to any MC 
simulation; however, in practice it is often the case that some sort of detector 
simulation is used to determine, for example, track-finding efficiencies. If a 
detector response were derived from data alone then this method would be 
completely data-driven and independent of any simulation. 

Another limit on the method is the quality of the fits to the iterations 
curves and the level of confidence in those fits that can be obtained. The 
fit encapsulates the behaviour of the detector; a bad fit to the iterations 
points will result in a bad description of the detector and a consequent poor 
correction. For this note we used a heuristic approach to finding the fitting 
function, which recognises that no function is a-priori more correct than any 
other; we simply wanted a function that described the behaviour of the de- 
tector. The performance of the fitted function can be estimated either by 
obtaining the original un-corrected data distributions from the additional 
applications of the detector effect or by performing a closure test on a sim- 
ulated sample of events; if the fit is good then the distributions obtained by 
passing hadron-level Monte Carlo events through a detector simulation and 
applying the correction procedure should agree with the original hadron-level 
results. Testing the procedure on a given observable and experimental setup 
by performing such closure tests on an ensemble of different Monte Carlo 
event samples is, in general, a good test of the method. 

A good correction also requires a reasonable description of the detector. 
The variation shown here between the nominal and alternative detector pa- 
rameterisations shows that the method has some degree of robustness against 
an imperfect parameterisation of the detector effect; however, as demon- 
strated by the effect of splitting, missing out an effect entirely can cause (or 
add to) a bias in the final distributions. 

The two particle correlation observable used here presents peaks in the 
distribution that are fixed in their Arj and A(p locations, regardless of the de- 
tector effect applied. Other observables may, however, exhibit peaks whose 
position is dependent on the detector effect (particle multiplicity being a sim- 
ple example). If such moving peaks are present in the observable then fitting 
an iterations curve independently to each bin of the observable may very 
well not be the optimum corrections strategy because it ignores information 
about the overall change in the shape of the distribution. In such cases it 
may be better to fit a function to the entire distribution and extrapolate the 
fit parameters to an A = — 1 detector effect in order to obtain a correction 
factor. We suggest that Bezier curves or their basis-spline generalisations 
may be ideal for this purpose; each additional application of the detector 
effect will cause the control points of the curve to migrate along a path. 

While the example presented here shows a charged-track based measure- 
ment, there is nothing that would prevent this method being extended to 
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measurements using calorimeter energy deposits. In order to perform a cor- 
rection on such an observable, a reasonable parameterisation of the energy 
and angular smearing produced in a calorimeter would be required. Exist- 
ing fast simulation packages may provide a suitable starting point for such 
parameterisations. Correcting calorimeter observables in this way may even 
allow jet energies to be calibrated and corrected for soft phenomena such as 
multi-parton interactions. 
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Figure 5: At] (a) and A<p (b) dependence of the two particle correlation func- 
tion corrected for the effects of particle splitting. Corrections were applied 
using polynomial fits of degree four to the iterations curves (see figure [3]) . 
The nominal detector-level observable (without splitting) is shown as filled 
circles, while the additional effect of splitting is shown as un-filled diamonds. 
Correcting the sample with splitting by using only the nominal detector pa- 
rameterisation (un-filled triangles) can lead to a significant bias. The effect 
of the splitting was completely removed (filled triangles) by performing the 
correction using a detector parameterisation that included the splitting ef- 
fect. The true value of the observable is shown as un-filled squares and the 
differences between the true and corrected values are shown in the bottom 
panel. 
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Figure 6: The uncorrected (iV = 0) two particle correlation function as 
observed by the nominal parameterised detector compared to values of the 
uncorrected function that were estimated only from fits to further (N > 1) 
applications of the detector effect. Fits to a polynomial of degree 3 (triangles) 
and degree 4 (circles) are shown for the Arj (a) and Acp (b) dependence of 
the two particle correlation function. 
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Figure 7: Polynomial fits to the loga- 
rithm of the absolute value of the two 
particle correlation function Vs. the 
number of applications of the detector, 
N. The value measured by the detec- 
tor is at N = 0, the true value is shown 
by a star at N = —1. The observed 
N = value is estimated from both a 
cubic and quartic fit to the N = 1 and 
higher points, while the true N = —1 
value is estimated from a quartic fit to 
the N = and higher points. Fits are 
shown for three A77 values of (a), 1.6 
(b) and 4.5 (c). 
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Figure 8: Corrected and uncor- 
rected A?y dependence of the two 
particle correlation function. Cor- 
rections were applied using poly- 
nomial fits of degree four to the 
logarithm of the iterations curves 
(see figure [7| . The Arj curves are 
otherwise as described in figure |4j 
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